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Speckle-turbulence  interaction  has  the  potential  for  allowing  single 
ended  remote  sensing  of  the  path  averaged  vector  crosswind  in  a  plane 
perpendicular  to  the  line  of  sight  to  a  target.  If  a  laser  transmitter  is 
used  to  illuminate  a  target,  the  resultant  speckle  field  generated  by  the 
target  is  randomly  perturbed  by  the  atmospheric  turbulence  as  it  propagates 
back  to  the  location  of  the  transmitter-receiver.  When  a  crosswind  is 
present,  this  scintillation  pattern  will  move  with  time  across  the  receiver- 
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20.  Abstract  (continued) 

/  A 

'aperture;  and  consequently,  the  time  delayed  statistics  of  the  speckle  field 
at  the  receiver  are  dependent  on  the  crosswind  velocity. 

A  continuous  wave  t'cw^/ laser  transmitter  of  modest  power  level  (a  watt 
or  two)  in  conjunction  with  optical  heterodyne  detection  has  been  used  to 
exploit  the  speckle-turbulence  interaction  and  measure  the  crosswind.  The  use 
of  a  cw  transmitter  at  10.6  microns  and  optical  heterodyne  detection  has  many 
advantages  over  direct  detection  and  a  double  pulsed  source  in  the  visible  or 
near  infrared.  These  advantages  include  the  availability  of  compact,  reliable 
and  inexpensive  transmitters;  better  penetration  of  smoke,  dust  and  fog; 
stable  output  power;  low  beam  pointing  jitter;  and  considerably  reduced 
complexity  in  the  receiver  electronics.  In  addition,  with  a  cw  transmitter, 
options  exist  for  processing  the  received  signals  for  the  crosswind  that  do 
not  require  a  knowledge  of  the  strength  of  turbulence. 

V\ 

During  the  period  of  this  grant,  an  optical  heterodyne  transmitter- 
receiver  system,  using  a  CO2  waveguide  laser  as  a  source,  was  designed  and 
built.  Several  important  results  were  obtained  with  the  system  including 
experimental  proof  that  for  coherent  detection  and  a  diffuse  target,  the 
peak  of  the  time  delayed  covariance  function  does  shift  with  crosswind  as  , 
required  for  crosswind  measurement.  Reciprocity  relationships  for  use  with  ? 
time  delayed  statistics  were  developed  and  verified,  and  the  crosswind  and 
strength  of  turbulence  was  remotely  sensed. 

Wind  and  turbulence  strength  measurements  were  made  at  500  meter  and 
1000  meter  ranges.  The  results  are  encouraging  and  have  very  positive 
implications  with  respect  to  remote  sensing  of  winds  and  turbulence  in  the 
atmosphere. 
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SUMMARY 


Speckle-turbulence  interaction  has  the  potential  for  allowing  single 
ended  remote  sensing  of  the  path  averaged  vector  crosswind  in  a  plane 
perpendicular  to  the  line  of  sight  to  a  target.  If  a  laser  transmitter  is 
used  to  illuminate  a  target,  the  resultant  speckle  field  generated  by  the 
target  is  randomly  perturbed  by  the  atmospheric  turbulence  as  it  propagates 
back  to  the  location  of  the  transmitter-receiver.  When  a  crosswind  is 
present,  this  scintillation  pattern  will  move  with  time  across  the  receiver 
aperture;  and  consequently,  the  time  delayed  statistics  of  the  speckle  field 
at  the  receiver  are  dependent  on  the  crosswind  velocity. 

A  continuous  wave  (cw)  laser  transmitter  of  modest  power  level  (a  watt  or 
two)  in  conjunction  with  optical  heterodyne  detection  has  been  used  to  exploit 
the  speckle-turbulence  interaction  and  measure  the  crosswind.  The  use  of  a  cw 
transmitter  at  10.6  microns  and  optical  heterodyne  detection  has  many 
advantages  over  direct  detection  and  a  double  pulsed  source  in  the  visible  or 
near  infrared.  These  advantages  include  the  availability  of  compact,  reliable 
and  inexpensive  transmitters;  better  penetration  of  smoke,  dust  and  fog; 
stable  output  power;  low  beam  pointing  jitter;  and  considerably  reduced 
complexity  in  the  receiver  electronics.  In  addition,  with  a  cw  transmitter, 
options  exist  for  processing  the  received  signals  for  the  crosswind  that  do 
not  require  a  knowledge  of  the  strength  of  turbulence. 

During  the  period  of  this  grant,  an  optical  heterodyne 
transmitter-receiver  system,  using  a  CO2  waveguide  laser  as  a  source,  was 
designed  and  built.  Several  important  results  were  obtained  with  the  system 
including  experimental  proof  that  for  coherent  detection  and  a  diffuse  target, 
the  peak  of  the  time  delayed  covariance  function  does  shift  with  crosswind  as 


required  for  crosswind  measurement.  Reciprocity  relationships  for  use  with 
time  delayed  statistics  were  developed  and  verified,  and  the  crosswind  and 
strength  of  turbulence  was  remotely  sensed. 

Wind  and  turbulence  strength  measurements  were  made  at  500  meter  and 
1000  meter  ranges.  The  results  are  encouraging  and  have  very  positive 
implications  with  respect  to  remote  sensing  of  winds  in  the  atmosphere. 
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Introduction 


Speckle-turbulence  interaction  has  the  potential  for  allowing  single 
ended,  remote  sensing  of  the  path  averaged,  vector  crosswind.  If  a  laser 
transmitter  is  used  to  illuminate  a  target  (Figure  1),  the  resultant  speckle 
field  gt  ated  by  the  target  is  randomly  perturbed  by  the  atmosphere  as  it 
propagates  back  to  the  receiver.  This  creates  a  scintillation  pattern  at  the 
receiver;  and,  if  a  crosswind  is  present,  this  pattern  will  move  across  the 
field  of  view  of  the  receiver.  Consequently,  the  time  delayed  statistics  of 
the  speckle  field  at  the  receiver  are  dependent  on  the  crosswind  velocity  and 
can  be  used  to  determine  the  vector  crosswind  in  a  plane  perpendicular  to  the 
propagation  direction. 

As  will  be  discussed  below,  we  have  previously  accomplished  some  work  on 
remote  crosswind  sensing  utilizing  speckle-turbulence  interaction,  a 
polychromatic  source  and  direct  detection. This  work  was  done  in  the 
near  infrared  utilizing  a  Nd:YAG  laser  operating  at  1.06  microns.  With  a 
direct  detection  receiver,  a  pulsed  transmitter  is  required;  and  furthermore, 
in  order  to  measure  the  time  delayed  statistics,  the  transmitter  must  be 
double  pulsed.  This  resulted  in  a  complicated  and  very  expensive  transmitter 
that  was  plagued  by  beam  alignment,  beam  jitter,  stability  and  reliability 
problems.  Nevertheless,  some  good  experimental  data  that  indicates  the 
speckle-turbulence  technique  has  merit  was  obtained. 

Examples  of  this  data  are  shown  in  Figures  2  and  3.  The  data  was  taken 
at  the  Biggs  Optical  Test  Range  (BOTR)  in  El  Paso,  Texas.  In  situ 
instrumentation  consisted  of  a  linear  array  of  20  propeller  anemometers  spaced 
evenly  along  the  path.  Their  outputs  were  added  without  weighting  to  provide 
the  solid  trace.  The  output  of  the  pulsed  system  is  shown  by  the  dashed 
lines.  A  four  second  averaging  time  was  used. 
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The  key  features  of  the  direct  detection  system  were: 

•  Pulsed  Laser 

•  Polychromatic 

•  Direct  Detection 

•  1.06  Microns  Wavelength 

Although  good  results  were  obtained  with  this  system,  it  has  a  fatal  flaw. 
Because  of  the  limited  pulse  repetition  rate,  the  only  usable  processing 
technique  is  the  slope  of  the  time  lagged  covariance  function  at  zero  time 
delay.  However  the  slope  is  proportional  to  both  the  crosswind  and  the 
strength  of  turbulence.  Since  the  strength  of  turbulence  is  not  known  a 
priori,  the  system  will  not  be  useful  until  there  are  significant  improvements 
in  pulsed  laser  technology  that  will  allow  other  signal  processing  methods  to 
be  used  that  do  not  require  an  a  priori  knowledge  of  the  strength  of 
turbulence . 

A  continuous  wave  (cw)  laser  transmitter  of  modest  power  level  in 
conjunction  with  optical  heterodyne  detection  can  be  used  to  exploit  the 
speckle-turbulence  interaction  and  measure  the  crosswind.  The  key  features  of 
such  systems  are: 

•  Continuous  Wave 

•  Monochromatic 

•  Heterodyne  Detection 

•  10.6  Microns  Wavelength 

which  is  just  about  the  opposite  of  the  features  for  the  1.06  micron  system. 

The  use  of  a  cw  transmitter  at  10.6  microns  and  optical  heterodyne 
detection  has  many  advantages  over  direct  detection  and  a  double  pulsed  source 
in  the  visible  or  near  infrared.  These  advantages  include  the  availabi  lity  of 
compact,  reliable  and  inexpensive  transmitters;  better  penetration  of  smoke, 
dust  and  fog;  stable  output  power;  low  beam  pointing  jitter;  and  considerably 


reduced  complexity  in  the  receiver  electronics.  In  addition,  with  a  cw 
transmitter,  options  exist  for  processing  the  received  signals  for  the 
crosswind  that  do  not  require  a  knowledge  of  the  strength  of  turbulence,  the 
fatal  flaw  of  the  1.06  microns  system.  During  the  period  of  this  research 
grant,  an  optical  heterodyne  transmitter-receiver  system,  using  a  C02 
waveguide  laser  as  a  source,  was  designed  and  built.  Several  important 
results  were  obtained  with  the  system  including  experimental  proof  that  for 
coherent  detection  and  a  diffuse  target,  the  peak  of  the  time  delayed 
covariance  function  does  shift  with  crosswind  as  required  for  crosswind 
measurement.  A  reciprocity  relationship  for  use  with  time  delayed  statistics 
was  developed  and  verified,  and  the  crosswind  and  strength  of  turbulence 
remotely  sensed. 

Wind  and  turbulence  strength  measurements  were  made  at  500  meter  and  1000 
meter  ranges.  The  results  are  encouraging  and  have  very  positive  implications 
with  respect  to  remote  sensing  of  winds  in  the  atmosphere. 

Background 

Analytic  formulations  for  the  first  and  second  order  statistics  of  the 
received  intensity  have  been  developed  for  the  case  where  a  diffuse  target  is 
illuminated  by  a  TEty0  monochromatic  laser  beam.9-H  They  properly  account 
for  the  effect  of  atmospheric  turbulence  both  for  propagation  of  the  laser 
beam  from  the  transmitter  to  the  target  and  propagation  of  the  target  induced 
speckle  field  back  to  the  receiver.  The  exact  formulations  are  quite  complex 
and  difficult  to  interpret.  However,  for  longer  wavelength  sources  such  as 
10.6  microns,  simple  formulations  based  on  the  joint  Gaussian  assumption  are 
valid  for  almost  all  physically  realizable  scenarios. 

From  previous  work,  the  time  delayed  covariance  function  using  the  joint 
Gaussian  assumption  is  given  for  the  focused  case  by? 
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(1) 
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2a  2  3p  5/3  1  1 
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where  P^  designates  the  location  of  a  detector  in  the  receiver  plane, 

P  =  P2  -  Pj 

T  =  t2  ~  ti 

o  =  Transmitter  Beam  Radius 
o 

=  Transverse  Phase  Coherence  Length 

V  =  Vector  Wind  Velocity 
and 

W  =  Normalized  Path  Length  from  the  Transmitter-Receiver  to  the  Target. 

As  can  be  seen  from  Eq.(l),  if  two  detectors  are  separated  by  a  vector 
distance  P,  the  time  delayed  covariance  can  be  measured;  and  from  Eq.(2)  it 
can  be  seen  that  the  measured  quantity  will  be  a  function  of  the  crosswind 
velocity . 

With  a  cw  system  there  are  several  options  for  processing  the  data  to 
obtain  the  wind.  Some  of  these  methods  are  illustrated  in  Figure  4.  The 
Briggs  method  measures  the  time  delay,  Tg,  at  which  the  autocovariance  and 
the  time-lagged  covariance  curves  cross.  The  delay  to  peak  method,  measures 
the  time  delay,  Tp,  where  the  time-lagged  covariance  reaches  its  peak 
value.  The  width  of  the  autocovariance  method  measures  the  time  delay,  xw, 
at  which  the  autocovariance  curve  decreases  to  67%  of  its  peak  value,  the 
slope  method  measures  the  slope,  S,  of  the  time-lagged  covariance  function  at 
zero  time  delay.  All  of  these  methods  have  been  used  for  the  line  of  sight 
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case  and  appear  to  have  applicability  to  remote  wind  sensing  using  speckle 
turbulence  interaction.  Assuming  uniform  wind  and  turbulence,  the  wind 
velocity  for  these  methods  is  given  by 


Briggs  Method 


V  = 


3.1056  t 


Delay  to  Peak  Method 


V  = 


2t  . 


Width  of  Autocovariance  Method  (67%) 


V_  = 
P_ 


0.1395 


Slope 


V  - 


3  po 


5  /3 


32  p2/3 


Cj(p,o) 


The  above  techniques  were  not  evaluated  during  the  period  of  this  report 
due  to  time  constraints.  However,  during  the  next  year  or  so,  it  is  planned 
to  evaluate  the  usefulness  of  each  of  these  techniques.  Some  new  methods  have 
been  developed  that  use  both  the  autocovariance  function  and  the  time  lagged 
covariance.  Formulation  for  these  quantities  are  functions  of  both  the 
crosswind  and  pD  the  transverse  phase  coherence  length.  By  combining  them 
and  also  averaging  the  results  for  several  different  time  delays,  both  the 
crosswind  and  pQ  can  be  measured.  One  result  for  the  crosswind  which  was 
used  to  process  the  data  presented  in  a  later  section  is 


(3) 
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where 

P  =  detector  spacing 

a  -  Beam  Radius 

o 

t ^  =  i**1  time  delay 

C_.,  *  Normalized  (to  the  mean  squared)  time  lagged  covariance 

IN 

N  =  Number  of  time  delays  used 

and  where  the  direction  of  the  crosswind  is  determined  by  the  skewness  of  the 
TLC.  There  are  other  variations  on  this  approach;  but  better  data  is  needed 
for  evaluation.  The  number  of  time  delays  used  for  the  data  processed  was 
typically  five. 


Experimental  Apparatus 

The  transmitter  optics  is  shown  in  Figure  5.  The  source  is  a  C02 
waveguide  laser.  Part  of  the  laser  radiation  is  split  from  the  main  beam  and 
focused  onto  the  two  detectors  for  use  as  an  optical  local  oscillator.  The 
remaining  part  of  the  beam  is  directed  through  an  acoustooptic  modulator  (AO) 
where  its  optical  frequency  is  shifted  by  37.5  MHz  and  then  passed  through  a 
AO  X  beam  expander  that  is  focused  on  the  target.  The  half-wave  plate  rotates 
the  polarization  of  the  beam  from  vertical  to  horizontal  to  match  the  AO 
modulator  and  the  quarter-wave  plate  provides  a  circularly  polarized  output 
beam. 

The  transmitter  beam  is  scattered  diffusely  by  a  4  foot  x  4  foot 
sandblasted  aluminum  target.  The  returning  radiation  is  directed  by  two 
one-inch  mirrors  onto  lenses  that  focus  it  on  the  detectors.  The  effective 
area  of  the  receivers  is  approximately  one-half  square  inch  and  the 


transmitted  power  is  around  one  watt. 
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The  main  problem  encountered  in  using  the  system  was  optical  feedback  of 
the  transmitted  beam  to  the  detectors.  The  source  of  the  feedback  is  the 
acoustooptic  modulator.  Apparently  there  is  some  diffuse  scattering  of  both 
the  unshifted  and  frequency  shifted  light  by  the  germanium  crystal  used  in  the 
modulator.  Even  though  there  is  around  100  db  of  isolation,  the  levels  fed 
back  are  sufficient  to  cause  large  homodyne  type  amplitude  beats  on  the 
heterodyne  signal.  The  feedback  light  causes  two  problems:  when  it  is 
detected  it  creates  a  DC  offset  that  makes  processing  for  the  statistics 
difficult;  and  it  also  mixes  with  the  actual  signal  from  the  target  which 
generates  large  amplitude  homodyne  beats. 

The  optical  feedback  problem  was  not  anticipated  and  a  great  deal  of 
effort  was  made  to  eliminate  it.  By  using  quarter-wave  plates  and  Brewster 
angle  polarizers  the  problem  was  reduced  to  a  point  wt  re  some  useful  but 
still  contaminated  data  could  be  obtained.  Unfortunately,  the  polarization  of 
the  optical  feedback  was  not  stable;  and  so  the  degree  of  isolation  would 
vary  over  time.  Frequent  readjustment  was  required. 

Two  solutions  were  proposed  for  the  problem:  a  gated  cw  transmitter;  and 
the  use  of  separate  lasers  for  the  transmitter  and  local  oscillator.  By  using 
the  time  delay  inherent  in  the  transmitted  wave  traveling  to  the  target  and 
back  to  the  receiver,  the  transmitter  and  receiver  can  be  synchronously  turned 
on  and  off  out  of  phase  such  that  the  optical  feedback  is  not  present  on  the 
receiver  at  the  same  time  as  the  signal  from  the  target.  The  laser  actually 
runs  cw  for  stability  reasons  and  the  transmitter  switching  is  done  with  the 
acoustooptic  modulator.  It  was  expected  that  two  to  three  orders  of  magnitude 
improvement  would  be  obtained.  However,  only  about  a  factor  of  25  was 
achieved  due  to  scattering  from  echoes  in  the  acoustooptic  (AO)  modulator. 

The  addition  of  a  second  AO  modulator  would  significantly  increase  the 


isolation . 


As  mentioned  above,  the  second  approach  uses  separate  lasers  for  the 
transmitter  and  LO.  Since  there  is  no  common  optical  path,  virtually  complete 
isolation  should  be  obtainable. 

The  method  of  gating  the  transmitter  and  receiver  is  illustrated  in 
Figure  6.  The  transmitter  was  gated  by  switching  the  RF  drive  to  the 
acoustooptic  modulator  and  the  receivers  were  gated  by  switching  the  RF  local 
oscillator  used  to  mix  down  the  received  signal  at  37.5  MHz  to  the  5  MHz  IF 
frequency.  A  block  diagram  of  the  receiver  is  shown  in  Figure  7.  The  data 
contained  in  the  next  section  was  taken  with  the  switched  cw  system  as 
described  above. 

Experimental  Results 

Representative  experimental  data  with  the  target  at  500  meters  range  is 
shown  in  Figures  8  through  10.  Equation  (3)  was  used  for  the  processing  and 
it  should  be  noted  that  there  are  no  arbitrary  parameters  used  and  that  no  a 
priori  knowledge  of  the  strength  of  turbulence  is  needed.  These  results  are 
encouraging . 

A  small  amount  of  data  has  been  obtained  with  the  target  at  1000  meters. 
Since  the  signal  is  reduced  at  longer  ranges;  but  the  optical  feedback  remains 
constant,  the  quality  of  the  data  decreases.  Figures  11  and  12  show  some  of 
the  better  sections  of  data  for  1000  meters  range.  It  is  noted  that  there  is 
some  bias  which  may  be  caused  by  residual  correlation  from  the  optical 
feedback. 

The  in  situ  data  for  both  the  500  meters  and  1000  meters  path  length  were 
obtained  by  using  a  Campbell  Scientific,  CA-9,  path  averaging  crosswind 
sensor.  It  should  be  noted  that  the  CA-9  and  Heterodyne  sensors  have 
different  path  weighting  functions  which  may  account  for  some  of  the 
differences  in  indicated  wind  speed. 
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Figure  12.  1000  meter  Experimental  Data 

Wind  =  5.0  m/s 
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ABSTRACT 

In  optical  propagation  problems  involving  atmospheric 
turbulence  and  single  ended  transmitter-receiver  systems, 
formulations  for  the  time-lagged  statistics  will  involve  point 
spread  functions  representing  optical  propagation  in  both 
directions.  In  order  to  use  these  formulations,  reciprocity  must 
be  used.  However,  for  time  lagged  statistics,  care  must  be 
exercised  when  using  reciprocity.  This  paper  shows  how  the 
reciprocity  principle  should  be  used  in  the  case  of  time  lagged 
statistics  and  presents  experimental  evidence  to  support  the 


Introduction 


In  single  ended  optical  propagation  problems  in  the  turbulent 
atmosphere,  the  transmitter  and  receiver  are  located  at  the  same 
end  of  the  propagation  path.  The  transmitted  radiation  propagates 
down  the  path  and  strikes  a  target  where  it  is  scattered  or 
reflected  back  towards  the  receiver.  This  radiation  interacts 
with  the  atmospheric  turbulence  on  both  the  down  and  back  paths. 
The  reciprocity  principle**^  says  that  if  a  point  source  at  r  in 
the  transmitter  plane  generates  a  certain  field  at  p  in  the 
receiver  plane;  and  if  the  point  source  is  moved  to  p,  it  will 
then  produce  that  same  field  at  r.  It  can  be  expressed  in  terms 
of  the  complex  wave  perturbation  due  to  a  point  source  in  the 
turbulent  atmosphere  as 

^(r,p,t)*4>'(p,7,t)  (1) 

where  the  unprimed  quantity  refers  to  propagation  from  the 
transmitter/receiver  plane  source  to  the  target  and  the  primed 
quantity  refers  to  propagation  from  the  target  plane  back  to  the 
transmitter/receiver  plane.  Clearly,  as  long  as  the  transit  time 
is  small  with  respect  to  the  decorrelat ion  time  of  the  atmosphere, 
the  reciprocity  principle  is  valid.  However  Eq.(l)  is  symbolic 
and  as  will  be  shown,  does  not  necessarily  imply  that  formulations 
for  ♦  and  can  be  used  directly  in  it  when  dealing  with  problems 
that  involve  atmospheric  propagation  in  both  directions  and  time 
lagged  statistics. 


N> ' 


The  complex  wave  perturbation  can  be  used  in  the  extended 
Huygens-Fresnel  formulation  to  determine  the  fields  at  the 
receiver,  including  the  effects  of  turbulence.  This  requires  that 
the  extended  Huygens-Fresnel  formulation  be  used  twice,  once  for 
the  path  down  and  once  for  the  path  back.^  Consequently,  the 
final  result  for  the  field  at  the  receiver  will  contain  both  ^  and 
\|> 1 .  When  this  expression  for  the  fields  is  used  to  evaluate  the 
statistics  at  the  receiver,  the  result  will  contain  statistics  for 
both  and  ’ . 

If  the  turbulence  is  isotropic,  and  homogeneous  in  the 
transverse  plane,  then  the  non-time  lagged  statistics  will  depend 
only  on  the  magnitude  of  the  vector  distance  between  field  points 
and  no  problems  arise  using  Eq.(l).^»^»^  However,  when 
time-delayed  statistics  are  considered  and  a  crosswind  is  present, 
they  will  in  general  depend  on  the  vector  distance  between  field 
points  and  the  statistics  become  anisotropic.  Since  vector 
results  sometimes  depend  on  the  direction  of  propagation  and  the 
coordinate  system,^  care  must  be  taken  in  using  reciprocity  when 
time-lagged  statistics  are  considered  and  propagation  in  two 
directions  are  involved. 

Figure  1  illustrates  an  example  of  the  problem  and  will  help 
lead  to  the  solution.  At  the  top  of  the  figure  is  shown  a 
correlation  type  of  crossvind  sensor  on  the  left  and  a  point 
source  on  the  right.  The  sensor  ha6  two  detectors  labeled  D2  and 
D1 .  The  time  lagged  correlation  then  will  peak  for  a  positive 
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value  of  t  *  Cj  -  and  the  instrument  will  indicate  the 
crossvind  velocity  as  V.  However,  if  the  crosswind  sensor  is  now 
rotated  about  the  y  axis  shown  in  Figure  2  and  exchanged  with  the 
point  source  as  shown  in  the  lower  part  of  Figure  1  it  will 
indicate  the  crosswind  as  -V  since  D2  now  sees  the  turbulence 
pattern  before  Dl.  To  straighten  things  out,  the  crosswind  sensor 
must  also  be  rotated  about  the  z  axis.  It  will  then  indicate  the 
correct  value  of  crosswind  but  the  sense  of  the  actual  wind  with 
respect  to  the  new  coordinate  system  will  be  reversed.  It  appears 
then  that  in  problems  of  this  type,  either  the  coordinates  need  to 
be  reversed  or  the  crossvind  reversed  in  order  for  \p  and  1  to  be 
interchangeable  through  reciprocity;  and  that  a  correct 
reciprocity  relationship  to  be  used  in  single  ended  atmospheric 
propagation  problems  involving  time  lagged  statistics  is 

<K7,  P,  t,  V)  *  ip  '  (p  ,  r,  t,  -V)  .  (2) 

At  this  point  some  readers  may  feel  that  Eq.(2)  is  obviously 
correct.  However,  for  the  skeptics,  an  experiment  has  been 
designed  that  has  clearly  different  predicted  results  depending  on 
whether  Eq.(l)  or  Eq.(2)  is  used.  First  it  will  be  shown 
analytically  in  the  next  section  that  Eq.(2)  is  correct,  and  then 
the  predicted  and  experimental  results  will  be  presented. 

However,  before  proceeding  it  should  be  emphasized  that  Eq.(2) 
does  not  deny  that  reciprocity  as  implied  symbolically  by  Eq.(l) 
is  not  valid.  It  only  states  that  when  the  two  propagation 
directions  are  mixed  in  a  problem,  that  the  wind  velocity  must  be 
reversed  in  the  reciprocity  relationship. 
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Analytical  Development 

Although  the  problem  being  considered  is  not  strictly  time 
harmonic,  the  time  harmonic  solution  to  Maxwell's  equations  that 
leads®  to  a  formulation  for  4<(r,  p,  t,  v)  can  still  be  used 
because  the  index  of  refraction  variations  due  to  the  atmospheric 
turbulence  are  small  and  vary  slowly  in  time-  The  formulation  is 
given  by 

<K~,  p",  t)=  exp[  ik-l~~i-  ] 


r  :J  Ir  -  r-  j2 


/  dxjdvjdtj  nj  (rj  ,  t)  exp[  ik  — ^ 


4  y  /[  21  (L  ~  21  ),J 


(3) 


where  nj  (rj  ,  t,  v)  is  the  random  part  of  the  index  of  refraction 
and  L  is  the  distance  from  the  transmitter  to  the  target.  If  the 
source  and  field  points  are  interchanged  in  (3)  and  it  is 
transferred  to  the  primed  coordinate  system,  it  becomes 


*  '(7,  p  ,  t)  -  — ■  exp[  ik-  22~  - 1  ] 


•  //  dxx  'dyj  'dtj  '  n \  '  (tx  ' ,  t)  exp[  ik(-i^~J - 


♦  lp  ~  n  ’.f.]  /[  (gl  *(l  -  zlt)]  u> 

2(L  -  *J  ') 


,  . 

.v  --.’ey  -%  .-n  ^ 


A  -V  A 


>V» 


It  can  be  seen  from  Eq.(4)  that  if  the  functional  form  for  ny '  is 
the  same  as  for  nj  ,  then  the  functional  form  for  \|«  will  be  the 
same  as  for  i|>  '  and  they  are  interchangeable.  Since  for  the 
problem  of  interest  ny  is  approximately  proportional  to  6T,  the 
temperature  deviation  from  the  mean,  the  differential  equation  for 
6T  will  be  used  to  determine  the  proper  form  of  reciprocity  to  be 
used  in  problems  involving  propagation  in  both  directions.  The 
differential  equation  for  6T  is  given  by9 


ill  ♦  v  2H  ♦  v  ill  ♦  v  ill 

3  t  x  3  x  y  5y  2  3  z 


XV26T 


(5) 


where  y  is  the  thermal  diffusivitv  and  V  =  Vx  +  Vy  +  Vz 

x  y  z 

represents  the  wind.  Transferring  Eq.(5)  to  the  primed  coordinate 
system  it  becomes 


3  6T 1  36  T '  3 6  T '  3  6 T ’ 

3 1  x  3  x  '  y  3  y '  z  3  z ' 


XV  ,26T' 


(6) 


It  can  be  seen  that  if  the  wind  field  is  reversed  in  Eq.(6), 
then  it  will  be  identical  to  Eq.(5)  and  consequently  the 
functional  form  for  nj  '  will  be  identical  to  that  for  ny  .  This  is 
equivalent  to  expressing  the  wind  field  in  the  primed  system.  The 
proper  relationship  to  be  used  then  appears  to  be 


*(7o,r ,t,V)  «  ♦,(7,ro,t,V’) 


«  <|i  ’(7,7o,t  ,-v) 


(7) 


In  the  next  section,  experimental  results  that  support  this 
contention  will  be  presented. 


Experimental  Results 

A  CC^  laser  transmitter,  a  diffuse  target  and  an  optical 
heterodyne  receiver  have  been  used  to  measure  the  time-lagged 
covariance  function  of  the  intensity  for  an  atmospheric 
propagation  path.  The  transmitter  and  receiver  are  located  on  one 
end  of  the  path  and  the  target  at  the  other  end.  Before 
presenting  the  experimental  results,  a  previously  derived 
formulation  for  the  time-lagged  covariance  will  be  used  to  show 
that  the  result  using  Eq . ( 1 )  is  clearly  different  from  the  result 
using  Eq . (2) . 

For  lower  levels  of  turbulence,  the  time-lagged  covariance 
for  the  focused  case  and  within  the  limitations  of  Taylor's 
hypothesis  is  given  by^ 

Cj  -  |r  (?!  ,P2  ,Vt  )|2  (e) 

where  T  is  the  mutual  intensity  given  by 


r(pl  ,P2  ,Vi  )  '(—j:)2  -y  exp^^2  -  P22)] 

it  Lr 

•  /  dTi  «*p[*p(pi  *  P)  "  7  ^.(P.o.Vt  )] 


JJ  drj  dr2  exp[  - 


r  (ri  +  r?  )  .  ik  —  — 

V - 1 — - — < —  4  r~  Pi  ‘  r 

2b  2 
o 


~2  D*  (o.T.V,) 


and  w'here  P1  and  P2  are  the  location  of  the  two  detectors  in  the 


transmitter /receiver  plane,  P  *  P2  -  Pj  ,  r  =  r2  -  rj  and  the  wave 
structure  functions  Ify  and  '  correspond  to  propagation  from 
the  transmitter  to  the  target  and  vice  versa  respectively.  The 
wave  structure  function  is  given  by 


D^(p  ,r  ,Vt  ) 


<exp[ 4/  (p2,r2,t2)  +  l^Cpj  ,rj  ,t!  )]  > 


2.91  Lk2  /  C  2  (w)  wp  +  (1-v)  7  - 
4  n 


T  _  77.  15/3 


In  developing  Eq.(9)  r,  p,  and  P  were  used  to  designate  the 
coordinates  in  the  transmitter,  target  and  receiver  planes 
respectively  and  were  all  expressed  in  the  transmitter  coordinate 
system.  It  should  be  noted  that  the  coordinates  at  the  target 
were  not  reversed  after  scattering  from  the  target  in  deriving 
Eq . ( 9) . 
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Eq.(9)  can  be  evaluated  by  first  performing  the  dpj  integration. 
This  yields  after  making  the  change  of  variables 
r2  ~r  *  r  and  r2  *  r;  *  2R 
V  2  . 

r(Pj  ,P2  ,Vt)  =  exp[±£  (Pi2-P22)  -  D^,(P,o,Vt)] 
v  L 

•  /  dR  /  dr  i  (P+r)  exp[  -  — — r-  —  -  -  -|D^(o,r,Vr)]  (10) 

2a 

o 

where  use  has  been  made  of 

/  depi  exp[^-  Pj‘  (P-7)]  =  2n  Jo(^P]|p^7|)  (11) 

and 

2 

/  Pi  J0(|  PI  [P+TpdP!  -  2rr  (i)  6  (P-7)  (12) 

The  dR  integrations  are  now  easily  made  and  the  final  result  is 
2  2 

n?!  ,?2  .Vt)  -  -p-  exp[-5T  (pi2  •  p22 > 

-|^((P,o,Vr)  -  \  D^(o,-P,Vt)  -  p2/4ao2]  (13) 

Equation  (13)  will  be  evaluated  using  Eq.(2)  which  is  the  proposed 
relationship  between  ♦  and  and  also  Eq.(l).  Using  these 
expressions,  the  relationships  for  the  wave  structure  functions 


become 


(14) 


D,  (o,-P,Vr  )  -  D,  ,(-P,o,-Vr  )  «  D,  ,(P,o,Vt) 

T  T  V 

and 

D^o.-P.Vr)  «  ^.(-P.o.Vr  )  -  D^,(P,o,-Vt)  (15) 

respectively.  The  tiire  lagged  covariance  is  plotted  in  Figure  3 
for  the  two  cases.  The  shape  of  the  curves  is  similar,  but  the 
curve  for  the  case  using  Eq.{2)  has  its  peak  shifted  from  the 
origin  in  the  direction  of  positive  t  for  positive  crosswind.  It 
would  shift  the  opposite  direction  for  negative  crosswind. 

However  the  curve  for  the  other  case  does  not  shift  with  the 
crosswind.  It  appears  then  that  shifting  or  not  shifting  of  the 
peak  of  the  time-lagged  covariance  with  the  crosswind  should 
provide  a  good  test  of  validity  for  the  use  of  Eq.(2)  in  single 
ended  time  lagged  problems. 

Figures  4  and  5  show  typical  experimental  data  for  the 
time-lagged  covariance  with  the  crosswind  in  opposite  directions. 
As  can  be  seen,  the  peaks  shift  as  predicted  for  the  case 
corresponding  to  Eq.(2)  and  this  along  with  the  analysis  leads  us 
to  conclude  that  Eq.(2)  is  the  proper  relation  to  be  used  in 
problems  of  this  type. 
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Figure  Captions 

1.  Wind  sensing  using  a  correlation  detector  in  the  line  of 
sight  case. 

2.  Path  geometry  and  the  coordinate  systems  used. 

3.  Normalized  time  delayed  covariance  of  the  intensity  of  a 
target  generated  speckle  pattern  versus  time  delay  in  the 
turbulent  atmosphere  with  and  without  wind  reversal;  Theory. 

4.  Normalized  time  delayed  covariance  of  intensity  of  a  target 
generated  speckle  pattern  versus  time  delay  in  the  turbulent 
atmosphere;  Experimental  data  using  a  CC^  heterodyne  system. 


5. 


Same  as  4. 
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Speckle-turbulence  interaction  has  the  potential  for  allowing  single 
ended  remote  sensing  of  the  path  averaged  vector  crosswind  in  a  plane 
perpendicular  to  the  line  of  sight  to  a  target.  If  a  laser  transmitter  is 
used  to  illuminate  a  target,  the  resultant  speckle  field  generated  by  the 
target  is  randomly  perturbed  by  the  atmospheric  turbulence  as  it  propagates 
back  to  the  location  of  the  transmitter-receiver.  When  a  crosswind  is  present, 
this  scintillation  pattern  will  move  with  time  across  the  receiver  aperture; 
and  consequently,  the  time  delayed  statistics  of  this  turbulence  perturbed 
speckle  field  are  dependent  on  the  crosswind  velocity.  A  continuous  wave  (cw) 
laser  transmitter  of  modest  power  level  in  conjunction  with  optical  heterodyne 
detection  can  be  used  to  exploit  the  speckle-turbulence  interaction  and  measure 
the  crosswind.  The  use  of  a  cw  transmitter  at  10.6  microns  and  optical 
heterodyne  detection  has  many  advantages  over  previous  work  which  used  direct 
detection  and  a  double  pulsed  source  in  the  near  infrared.  These  advantages 
include  the  availability  of  compact,  reliable  and  inexpensive  transmitters; 
better  penetration  of  smoke,  dust  and  fog;  stable  output  power;  no  beam 
pointing  jitter;  and  considerably  reduced  complexity  in  the  receiver  electronics. 
In  addition,  with  a  cw  transmitter,  options  exist  for  processing  the  received 
signals  for  the  crosswind  that  do  not  require  a  knowledge  of  the  strength  of 
turbulence.  Analytical  and  experimental  results  will  be  presented. 

I-  I  ^This  material  is  based  in  part  upon  work  supported  by  the  Army 

I  Research  Office  under  contract  DAAG29-83-K-0077. 
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Abstract 

Atmospheric  turbulence  and  target  induced  speckle  can  have  a  significant  and  deleterious 
effect  on  the  performance  of  optical  systems.  The  primary  effect  is  to  introduce  a  strcr.r 
fluctuation  in  the  intensity  and  a  random  perturbation  of  the  phase  of  the  fields  at  the 
optical  receiver.  This  is  equivalent  to  introducing  a  large  noise  source  at  the  receiver; 
and  consequently  multiple  samples  and/or  clever  processing  techniques  are  needed  to  accu¬ 
rately  estimate  the  magnitude  of  the  received  intensity.  In  addition,  for  coherent  detec¬ 
tion  schemes,  the  random  phase  and  finite  transverse  coherence  length  introduced  by  the 
speckle  and  turbulence  will  limit  the  size  of  aperture  that  car.  be  used.  These  character¬ 
istics  of  the  speckle  field  after  propagation  through  turbulence  and  its  effect  on  optical 
systems  will  be  discussed  in  a  tutorial  overview. 


Z  ntrod ucti on 


Atmospheric  turbulence  and  target  induced  speckle  can  have  a  significant  and  deleterious 
effect  on  the  performance  of  optical  systems  that  operate  in  or  through  the  turbulent 
atmosphere.  Figure  1  shows  a  typical  system  arrangement. 


Figure  1.  Typical  system  configuration. 


Figure  2.  Signal  contaminated  by 
speckle  with  unity  contrast  ratio. 


The  source  is  s  TEM^  Zsser  beam  that  propsgstss  through  the  atmospheric  turbulence  to  n 
diffuse  target.  The  turbulence  will  degrade  the  spatial  coherence  of  the  beam,  in¬ 

troduce  amplitude  fluctuations  and  cause  the  beam  to  wander  over  the  target  surface.  The 
diffuse  surface  will  sonpletly  destroy  the  spatial  coherence  of  the  beam  and  generate  a 
speckle  field  that  will  then  propagate  back  through  the  turbulence  to  the  receiver.  The 
primary  effects  at  the  receiver  are  a  strong  fluctuation  in  the  intensity  and  a  random, 
perturbation  of  the  phase  of  the  fields.  This  is  squivalent  to  introducing  a  large  noise 
source  at  the  receiver  end  consequently  multiple  samples  and/or  clevsr  processing  tech¬ 
niques  are  needed  to  accurately  estimate  the  magnitude  of  the  received  intensity.  In 
addition,  for  coherent  detection  schemes,  the  random  phase  and  finite  transverse  coherence 
length  Introduced  by  the  speckle  and  turbulence  will  limit  the  size  of  aperture  that  can  be 
meed. 


Background 


Before  discussing  the  combined  effects  of  speckle  end  turbulence,  some  background  on 
speckle  atatistics  will  be  presented  to  give  the  reader  e  better  feel  for  the  speckle  prob¬ 
lem  before  the  combined  effects  of  speckle  and  turbulsnee  are  introduced.  The  statistical 
parameters  that  are  isportant  ere  the  normalised  variance,  which  ia  e  measure  of  the  degree 
of  signal  fluctuation,  the  covariance,  from  which  the  transverse  correlation  length  can  be 
derived,  and  the  probability  density  function.  For  a  monochromatic  laser  source,  a  diffuse 
target,  and  no  turbulsnee  the  normalised  variance  of  the  received  intensity  is  equal  to 
unity1  and  is  dafinsd  as 
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A  sample  function  for  chic  type  of  signal  is  shown  in  Figure  2  where  the  horizontal  axis  is 
time  or  space  and  the  vertical  axis  is  the  received  intensity.  Zn  the  former  case  there 
must  be  relative  motion  between  the  incident  beam  and  the  target.  Otherwise  the  signal  for 
any  particular  sample  function  would  be  constant  in  time.  There  would  however  be  a  sub¬ 
stantial  variation  in  signal  level  between  the  sample  functions  making  up  a  complete 
ensemble.  Equation  (1)  is  the  result  of  taking  an  ensemble  average.  However,  the  sane 
result  would  be  obtained  if  the  sample  function  in  Figure  2  were  averaged  to  obtain  the 
mean  and  variance.  In  any  case,  the  speckle  with  its  unity  normalized  variance  constitutes 
a  strong  source  of  random  noise  in  the  system. 


The  normalized  covariance  of  the  received  intensity  is  defined  as 


C.  (P) 
■‘n 


« (1 (P: ) 


0 


1>) (I (p. ) 

2 

I 


<!>)> 


(2) 


where  p  «  p;  -  p, 

and  in  the  absence  of  turbulence  can  be  expressed  as2 

V?  -  .*p  j-  V  [77*  [t  -.0  -  *)]]j 


where  a  is  the  illuminating  beam  radius  at  the  transmitter,  p  is  the  spacing  between  de¬ 
tectors?  L  is  the  path  length,  F  is  the  distance  at  which  the  transmitter  is  focused  and  k 
is  the  wave  number.  It  can  be  seen  that  for  a  focused  transmitter,  the  speckle  size  back  at 
the  receiver  is  approximately  the  same  size  as  the  transmitted  beam.  As  the  focus  distance 
F  is  increased,  the  spot  size  on  the  target  grows  and  the  speckle  size  back  at  the  receiver 
decreases  continuously. 

The  speckle  scale  aize  is  important  if  aperture  averaging  is  used  to  reduce  the  variance 
of  the  received  signal.  Each  speckle  can  be  thought  of  as  being  independent  of  all  the 
other  speckles.  Consequently,  what  the  collection  aperture  does  then  is  add  up  the  inten¬ 
sities  from  each  speckle  and  since  they  are  independent,  the  variance  goes  down  as  1/t; 
where  N  is  the  number  of  independent  speckles.  This  effect  can  be  estimated  using  the  co- 
variance  scale  size  and  the  aperture  size  to  find  an  approximate  N;  or  an  exact  calcula¬ 
tion  can  be  performed. 1 


Turbulence  and  Speckle 

When  turbulence  is  present,  an  additional  fluctuation  is  introduced  in  the  signal  at  the 
receiver  in  addition  to  that  caused  by  speckle.  This  is  illustrated2  in  Figure  3  which 
shows  the  normalized  variance  of  the  received  intensity  versus  the  Rytov  variance  which  is 
the  normalised  variance  as  observed  by  a  point  detector  located  at  the  receiver  due  to  a 
point  source  at  the  target.  As  can  be  seen  from  the  figure,  the  normalized  variance  in¬ 
creases  up  to  a  Rytov  variance  of  around  0.1  and  then  decreases  until  it  reaches  unity 
again.  The  maximum  increase  is  around  25  percent.  Detailed  formulations  from  which  the 
variance  can  be  calculated  are  contained  in  Reference  2.  It  should  be  noted  however  that 
the  results  for  different  path  lengths,  beam  sizes  and  wavelengths  do  not  differ  much  from 
that  ahown  in  Figure  3.  The  Rytov  variance  is  given  by 

e  *  •  0.124  kV*  L11/*  C  J  <4) 

x  n 

where  CB*  represents  the  strength  of  turbulence  (structure  parameter) . 

There  ie  an  additional  effect  that  can  increase  the  variance  beyond  that  predicted  in 
Figure  3.  If  the  beam  wanders  on  and  off  the  target  or  wanders  around  on  a  target  of 
greatly  varying  reflectivity,  the  variance  can  increase  sisrkedly. *  This  wander  can  be 
caused  by  the  atsospheric  turbulence  or  by  beam  jitter  in  the  transmitter  system. 

In  order  to  evaluate  the  effects  of  aperture  averaging,  the  decorrelation  due  to 
separation  of  staples  In  apace  and/or  time,  and  the  transverse  phase  coherence  length, 
the  tine  delayed  eovarlance  is  needed.  At  lower  values  of  integrated  path  turbulence 
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Figure  3.  Normalized  variance  of  the 
received  intensity  versus  the  log- 
amplitude  variance,  focused  beam,  500-m 
range,  1.06-  and  10. 6-. m  wavelengths, 
a  *  3.795  cm. 
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Figure  4.  Normalized  covariance  of  the 

received  intensity,  focused  beam,  2000-n 

range,  10. 6-. m  wavelength,  a  =  3.795  cr. 
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(exJ  <  0.003),  a  very  aimple  analytic  formulation  for  the  time  delayed  covariance  can  be 
used.  It  is  given  for  the  case  of  uniform  turbulence  along  the  path  by  *• 
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where  c  the  transverse  phase  coherence  length  is  given  by 


o  *  (0.545625  C  2  Lk2)”3/5 
o  n 


V  is  the  vector  crosswind,  and  t  is  the  time  delay  between  samples.  It  should  be  noted 
that  Equation  (5)  does  not  include  the  effects  of  the  log-amplitude  perturbation  and 
therefore  is  not  complete.  Consequently  some  care  should  be  exercised  in  its  use.  For 
Instance,  Equation  (5)  predicts  that  the  normalized  variance  (C.  (0,0)/<I>2)  is  unity 
independent  of  turbulence  level  whereas  the  more  complete  formulation2  which  includes  the 
log-amplitude  perturbation  gives  the  result  (verified  by  experiment)  shown  in  Figure  3.  A 
sore  complete  but  considerably  more  complicated  formulation  is  available2  that  includes  the 
log-amplitude  perturbation.  Using  it,  numerical  data  for  the  covariance  curves  shown  in 
Figure  4  was  generated.  As  discussed  in  the  previous  section,  there  is  a  target  induced 
speckle  acale  size  that  for  a  given  scenario  is  fixed  CIO” 1  *  curve).  However,  when  turbu¬ 
lence  is  present,  this  scale  size  is  dynamic  and  as  illustrated  in  Figure  4  decreases  as 
the  strength  of  turbulence  increases. 

Figure  5  shows  experimental  and  numerical  data  (using  the  complete  theory)  for  the  time 
delayed  covariance  function, *»*  which  can  be  used  to  estimate  the  degree  of  correlation 
between  settles  that  are  separated  in  space  and  time.  Using  Equation  (5)  (o  2  «  0.003), 
the  autocovariance  for  the  received  intensity  can  be  obtained.  It  is  given  oy 


Cj(o,Vt)  ■  <I>2  exp 


-10.(7 


eri 


It  should  be  noted  that  both  the  time-delayed  covariance  and  the  autocovariance  depend  on 
the  crosswind  velocity  and  the  strength  of  turbulence.  The  Fourier  transform  of  Equation 
(()  yields  the  power  spectral  density  for  the  received  intensity  which  is  plotted  in  Figure 


Figure  5.  Time  delayed  covariance,  Figure  6.  Temporal,  power  spectral  density,  received 
numerical  and  experimental  data.  intensity. 


6.  It  can  be  used  to  estimate  the  bandwidth  requirements  for  experimentally  measuring  the 
speckle-turbulence  statistics  (0  2  <  0.003). 

In  estimating  the  effect  that  turbulence  and  speckle  will  have  on  system  performance, 
the  probability  density  function  (PDF)  is  useful.  For  a  monochromatic  source,  a  diff-sa 
target,  and  no  turbulence,  at  the  receiver  the  fields  have  a  Gaussian  distribution,  the 
field  amplitude  has  a  Rayleigh  distribution,  and  the  intensity  has  an  exponential  distri¬ 
bution.  Turbulence  will  modify  these  distributions.  Ho/ever,  for  integrated  path  turbu¬ 
lence  corresponding  to  Rytov  variances  less  than  0.003,  they  are  a  reasonable  approxi¬ 
mation.  Figure  7  shows  experimental  cumulative  PDF  data  for  Rytov  variances  of  0.0022  and 
0.1735.  At  the  lower  value,  the  data  matches  the  experimental  distribution  quite  well. 
However  at  the  hiaher  turbulence  level  there  is  substantial  deviation  from  exponential. 
Analytical  and  experimental  work®,  9  has  shown  that  for  a  fully  developed  Gaussian  speckle 
Pattern,  the  PDF  of  the  intensity  changes  from  an  exponential  density  function  at  zero 
turbulence  to  a  K  distribution  as  the  turbulence  level  increases.  The  PDF  is  given  by 

M+l  M-l 

Pj(o)  ■  Vi  (2  Vt)  (7) 


where  is  the  modified  Bessel  function  of  order  M-l,  i  equals  «I>,  and  M  equals 

2/(Oj  3  -  1)  Figure  8  shows  numerical  data  using  Equation  (7)  and  experimental  data  for 

the  cumulative  PDF.  Good  agreement  was  obtained.  The  parameter  M  can  be  obtained  in  terms 
of  the  path  length,  wavelength,  beam  size  and  strength  of  turbulence  using  Reference  2. 
Formulations  for  the  PDF  with  turbulence  effects  are  also  available8  for  partially 
developed  speckle,  speckle  with  a  coherent  background  (glints),  and  polychromatic  speckle. 


Beam  jitter  can  contribute  both  to  tenoral  fluctuations  of  the  signal  and  to  decorrela¬ 
tion  between  signal  samples.  A  measure  of  the  magnitude  of  this  effect  can  be  obtained 
from  the  crosscorrelation,  C,(o,t),  between  samples  separated  in  time.  Figure  9  shows  some 
experimental  data  taken  using  a  pulsed  Nd: TAG  laser. 

Tbs  crosscorrelation  (normalised  to  <X>3)  was  measured  for  the  received  signals  from  two 
successive  pulses  spaced  one  millisecond  apart.  Xf  the  beams  corresponding  to  the  two 
pulses  ere  perfectly  aligned,  then  they  will  go  through  the  same  portion  of  the  atmosphere, 
strike  the  asms  spot  on  the  target  and  generate  s  unity  normalised  erosseovariance  back  at 
the  receiver.  However,  if  the  beams  have  a  static  or  dynamic  misalignment,  the  signals 
will  decorrelate  and  the  crosscovariance  will  decrease  to  a  value  less  than  unity.  As  can 
be  seen  from  Figure  9,  there  appears  to  be  considerable  static  and  dynamic  misalignment  of 
the  laser  transmitter  used  in  the  experiment.  This  was  confirmed  by  experiments  on  the 
laser  transmitter.  The  misalignment  problem  has  been  analysed  and  formulations  are  avail¬ 
able10  which  Include  the  effects  of  static  and  dynamic  beam  misalignment.  Figure  10  shows 
some  numerical  data  using  these  formulations  for  the  dynamic  misalignment  case.  The 
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parameter  e  ia  defined  to  be  the  ratio  of  the  beam  misalignment  angle  corresponding  to  84 
percent  cumulative  probability  to  the  beam  divergence  angle  in  vacuum. 


Figure  7.  Probability  distribution  functions  Figure  8.  Comparison  of  experimental 
of  the  received  irradiance  seen  in  two  turbu-  cumulative  PDF  with  theory  for  a  mono- 

lence  conditions.  (a)  0.  experimental  data  chromatic  speckle  pattern,  300-m  path 

from  low  turbulence:  C  2  »  1.0  *  10' ls  m*2'3,  length. 

L  «  500  m,  c_  2  •  1.02.  (b)  •,  experimental 

N 

data:  Cn2  m  8.0  *  10" 1 ^  m"2'3,  L  ■  500  m, 
o.  2  -  1.2.  The  solid  lines  indicate  expo- 

nential  distributions  of  the  same  means  as 
the  experimental  points. 


It  should  be  noted  that  changing  the  wavelength  between  smasurements  can  have  the  same 
effect  as  moving  the  spot  on  the  target  in  that  it  changes  the  speckle  pattern  and  de¬ 
correlates  the  signals.  This  is  a  target  dependent  phenomenon;  and  in  general  if 
tu  <<  e.,  where  o  is  the  standard  deviation  of  the  target  surface  roughness  and  Aw  is  the 
change  in  frequency,  then  the  pattern  does  not  change  very  much  and  the  signals  will  be 
highly  correlated. 
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Figure  9.  Crosscovariance,  Target  Range  “ 
1600  m.  Time  Delay  ■  10' 3  a,  Average 
Crosewind  ■  3  n/» ,  average  o  “  0.38,  i  » 
1.06  microns,  40  samples  perxdata  point. 
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Figure  10.  Normalized  autocovariance  of  a 
randomly  misaligned  laser  beam  versus  C  J. 
Focused  beam,  500-m  range,  1.06- urn  wave¬ 
length,  e  ■  3.795  cm. 
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In  single  ended  optical  propagation  problems,  the  transmitter  and  receiver 
are  located  at  the  same  end  of  the  propagation  path.  The  transmitted  radiation 
propagates  down  the  path  and  strikes  a  target  where  it  is  scattered  or  reflected 
back  towards  the  receiver.  The  radiation  interacts  with  the  atnrcspheric 
turbulence  on  both  the  down  and  back  paths.  In  analyzing  this  type  of  problem, 
a  reciprocity  relationship  for  the  point  spread  function  is  very  useful.  For 
the  non-time  lagged  case,  there  is  a  well  known  reciprocity  relationship  and 
it  would  seem  that  Taylor's  Hypothesis  could  be  used  to  convert  it  to  the  time 
lagged  case.  However,  for  the  non-time  lagged  case,  the  statistics  are 
isotropic;  but  for  the  time  lagged  case  the  statistics  depend  on  direction. 
Consequently,  careful  consideration  must  be  given  to  the  coordinate  system 
used  to  describe  the  propagation  on  each  path.  A  time  lagged  point  spread 
function  will  be  developed  and  experimental  evidence  presented  to  support 
the  result. 
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THE  EFFECTS  OF  TARGET  AND  SOURCE  MOTION  AND  BEAM 
SLEWING  ON  SPECKLE  PROPAGATION  THROUGH 
TURBULENCE1 
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Atmospheric  turbulence  and  target  induced  speckle  can 
have  a  significant  and  deleterious  effect  on  the  performance 
of  optical  systems.  One  method  that  is  used  to  overcome 
these  effects  is  to  take  averages  in  space  and/or  time. 
However,  the  reduction  in  signal  fluctuation  due  to  averaging 
are  affected  by  target,  receiver  and  source  motion  and  by 
beam  slewing.  Consequently,  formulations  for  the  mutual 
intensity,  the  time  delayed  covariance  and  the  power  spectral 
density  that  include  these  effects  are  needed.  Formulations 
for  these  statistical  measures  have  been  developed  that 
include  the  effects  of  source,  receiver  and  target  motion  and 
beam  slewing  as  well  as  the  effects  of  the  speckle  and 
turbulence.  Utilizing  the  extended  Huygens-Fresnel 
formulation,  expressions  for  the  mutual  intensity,  the  time 
delayed  covariance,  and  the  power  spectral  density  have  been 
developed  that  are  valid  at  all  levels  of  turbulence.  A 
TEM00  monochromatic  laser  in  conjunction  with  a  perfectly 
diffuse  target  is  used  to  generate  the  speckle  field.  The 
formulations  obtained  are  reduced  to  numerical  data  for 
illustrating  the  magnitude  of  the  effects. 
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Atmospheric  turbulence  and  target  induced  speckle  can  have  a  significant 
and  deleterious  effect  on  the  performance  of  optical  systems.  One  method  that 
is  used  to  overcome  these  effects  is  to  take  averages  in  space  and/or  tire. 
However,  the  reduction  in  signal  fluctuation  due  to  averaging  are  affected 
by  target  and  source  motion.  Consequently,  formulations  for  the  mutual  inter.sit 
the  time  delayed  covariance  and  the  power  spectral  density  that  include  target 
and  source  motion  effects  are  needed.  Formulations  for  these  statistical 
measures  have  been  developed  that  Include  the  effects  of  source  and  target 
motion  as  well  as  the  effects  of  the  speckle  and  turbulence.  Utilizing  the 
extended  Huygens-Fresnel  formulation  and  assuming  that  the  optical  fields  are 
jointly  Gaussian  at  the  receiver,  expressions  for  the  mutual  intensity,  the 
time  delayed  covariance,  and  the  power  spectral  density  have  been  developed. 

A  TEMqo  monochromatic  laser  in  conjunction  with  a  perfectly  diffuse  target  is 
used  to  generate  the  speckle  field.  The  formulations  obtained  are  easily 
reduced  to  numerical  data,  but  are  valid  only  for  Rytov  variances  less  than 
0.03  due  to  the  jointly  Gaussian  assumption. 
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